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Consider an ensemble of N x N non-Hermitian matrices in which 
all entries are independent identically distributed complex random 
variables of mean zero and absolute mean-square one. If the entry 
distributions also possess bounded densities and finite (4 + e) mo- 
ments, then Z. D. Bai [Ann. Probab. 25 (1997) 494-529] has shown 
the ensemble to satisfy the circular law: after scaling by a factor of 
1/y/N and letting N — > oo, the empirical measure of the eigenval- 
ues converges weakly to the uniform measure on the unit disk in 
the complex plane. In this note, we investigate fluctuations from the 
circular law in a more restrictive class of non-Hermitian matrices 
for which higher moments of the entries obey a growth condition. 
The main result is a central limit theorem for linear statistics of 
type Xn{/) = X/feLi wnere Ai, A2, . . . , Ajv denote the ensem- 

ble eigenvalues and the test function / is analytic on an appropriate 
domain. The proof is inspired by Bai and Silverstein [Ann. Probab. 
32 (2004) 533-605], where the analogous result for random sample 
covariance matrices is established. 

1. Introduction. Eigenvalues of random Hermitian matrices in the limit 
of large dimension have been widely studied on account of their relevance to 
multivariate statistics and theoretical physics. The study of non-Hermitian 
random matrix ensembles remains undeveloped by comparison. While con- 
sidered as early as the mid-1960s by Ginibre [15], only in the last several 
years has the non-Hermitian case been picked up with serious interest, no- 
tably in the physics community ([10, 11] and [24] offer guides to that liter- 
ature). 

Regardless of symmetry considerations, the most basic object of study 
in any random matrix ensemble is the density of states, or the limiting 
spectral distribution of the eigenvalues. In the setting considered here, this 
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may be described as follows. Consider an N x N matrix in which all entries 
are independent identically distributed and standardized to have mean zero 
and mean-square one. If this matrix is normalized by a factor of 1/y/N 
and the limit N — > oo taken, then (modulo, perhaps, additional technical 
assumptions) it is anticipated that the counting measure of eigenvalues will 
converge weakly to the uniform measure on the unit disk in the complex 
plane C. This is referred to as the circular law. 

In the fundamental work [15], the cases in which all entries are either inde- 
pendent identically distributed complex or quaternion Gaussians are shown 
to be integrable. That is, the full joint density of the N eigenvalues is com- 
puted. Employing this explicit eigenvalue density, the first proper proof of 
the circular law for the complex Gaussian ensemble (Ginibre ensemble) may 
be found in [19]; see also [9] which contains a weaker form of the circular law 
for ensembles of real Gaussian entries. Using an important idea of Girko [14], 
Bai [1] has proved the most general result to date: the circular law holds 
for ensembles of independent complex entries such that the common entry 
distribution possesses a bounded density and finite (4 + e) moment. 

This paper is concerned with fluctuations about the circular law. We focus 
on linear spectral statistics of the form 

N 

(1-1) **(/) = £/(A fc ), 

k=l 

where Ai, A2, . . . , \n denote the ensemble eigenvalues. Note that Xj^(f) = 
N J /(A)/xjv(^A) where fijy(d\) is the empirical spectral distribution. For 
ensembles of Hermitian matrices, there is an enormous body of work on 
Gaussian distributional limits of X^{f) as N — ► 00, for which a variety of 
methods has been developed. We cite [3, 7, 16, 20, 23, 27] and [28] as a 
representative list. We should also mention a different class of ensembles for 
which similar investigations have been carried out. These are the classical 
groups: N x N unitary, orthogonal or symplectic matrices distributed ac- 
cording to Haar measure. While certain techniques may be imported from 
the Hermitian setting, these ensembles have their own structure and their 
analysis requires a separate set of tools; see, for example, [8, 18, 21, 22] 
and [30]. One common theme in all of these results is how they manifest 
the rigidity of the point processes formed by random matrix eigenvalues. 
For smooth test functions /, the linear statistic need only be centered to 
produce a central limit theorem, no normalization is required. [Even for dis- 
continuous test functions, the growth of the variance of Xjy(f) is much less 
than that for a set of N independent particles.] 

Thus far, fluctuation results for linear statistics of non-Hermitian random 
matrices have been relegated to Ginibre's complex ensemble, where the ex- 
plicit joint density of eigenvalues may be exploited. In [10], advantage is 
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taken of the particularly simple form of the density of the eigenvalue moduli 
in the Ginibre ensemble and a central limit theorem is described for statis- 
tics Xn(/) — E[Xpf(f)], where /(A) =/(|A|) and / is differentiable. [Again 
the variance is 0(1) as N — > oo.] Reference [26] extends this picture to dis- 
continuous statistics of the eigenvalue moduli and also proves a collection of 
results for the case where / depends only on the eigenvalue phase. Here, we 
are able to go beyond the Gaussian (Ginibre) setup and establish a central 
limit theorem for spectral statistics (1.1) in a class of non-Hermitian random 
matrices built out of more general entry distributions. The proof is inspired 
by the work of the second author and Z. D. Bai [3]. The method is fairly 
robust in terms of the underlying distribution of the matrix entries but does 
require the analyticity of the test function / in X^(f). 

The particular ensemble we consider, denoted by M, may be described as 
1/>J~N x a base N x N matrix M. The latter is comprised of independent 
identically distributed complex entries m« which are normalized in the usual 
way, 

£[mu]=0 and E[|mn| 2 ] = l 

and which also satisfy 

(i) E[m 2 u ]=0, 
(1.2) (ii) E[\rh u \ k ] < k ak for k > 2 and some a > 0, 

(hi) 3?(mn) and 3?(mn) possess a bounded joint density. 

Condition (1.2.i) is a Gaussian-like assumption; it is analogous to require- 
ments (ii) and (hi) in Theorem 1.1 of [3]. The moment condition (1.2.ii) 
supplies a priori control on the magnitude of the eigenvalues. In particular, 
(1.2.ii) implies that with probability one as N — > oo, the spectral radius of M 
converges to 1 and the spectral norm of M converges to 2, see [12] and [13]. 
While it is more natural to consider the spectral radius in the non-Hermitian 
setting, it will be the spectral norm which figures more prominently below. 
Finally, the smoothness condition (1.2. hi) is of a technical nature. The same 
assumption was put to similar use in the proof of the circular law in [1]. 

With the ensemble defined we may state our hrst result: 

Theorem 1.1. Let the independent entries of the non-Hermitian matrix 
M satisfy the conditions outlined in (1.2). Consider test functions f±, f%, . . . , 
analytic in a neighborhood of the disk \z\ < 4 and otherwise bounded. Then 
as N — > oo, the vector 



(1.3) (X N (h) - iY/i(0),Xjv(/ 2 ) - Nf 2 (0), ■ ■ .,X N (f k ) - Nf k (0)) 
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converges in distribution to the mean-zero multivariate Gaussian (X(/i),X(/2), 
X(fk)) with covariances 

(1.4) E[X(f e )X(U] = - [ ±f e ( z )±f m ( z )d 2 z, 

7r Ju dz dz 

in which U is the unit disk and d 2 z = d$l(z) d^s(z). 

Note that Xisr(f) has been centered about the asymptotic mean. For 
large N, E[X^{f)] is well approximated by iVx the average of / against 
the circular law and 

- [ f(z)d 2 z = f(0) 

if / is analytic. The limit distribution of X^(f) — E[X^{f)\ is the same, as 
will become apparent in the course of the proof. The assumed analyticity of 
the test functions is a by-product of the strategy of proof: our main accom- 
plishment in this note is actually a central limit theorem for the resolvent 
of the random matrix M. Central limit theorems for analytic statistics are 
then inferred by way of the Cauchy integral formula. More to the point, for 
any domain IcC on which / is analytic, 

(1.5) X N (f) - Nf(0) = / f(z){tr(z - M)" 1 - Nz' 1 } dz 

2m Jan 

holds on the event that all of Ai, . . . , Ajy are contained in D. Assuming that 
the probability of the complement of that event decays sufficiently fast 
[which is the content of estimate (2.1) below], Theorem 1.1 is a corollary 
of the next result. 

Theorem 1.2. View the centered resolvent 

<& N ( Z ) = tr(z - M)~ l - Nz' 1 

as a process in z G C, any fixed contour lying exterior to the disk of radius 
four. Then the family {<3n(z)} is tight in the space of continuous functions 
on C and converges weakly to the complex Gaussian process &(z) defined by 

(1.6) E[&{z)^w)] = {l- zwY 2 
and E[<5(z)\ =0. 

To pass the central limit theorem for &n(z) onto X^{f) requires only that 
we note the continuity of the map (5jy — > f c f(z)(5]y(z) dz from the space of 
continuous functions on C into C. As for connecting the covariance structures 
(1.4) and (1.6), choose functions / and g analytic in a neighborhood of 
U. Also, choose a contour C lying within the region of analyticity of both 
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functions, but enclosing U. Then (1.4), applied to X(f) and X(g), may be 
rewritten as 



(1.7) 



1 

7T Ju 



f'(v)g'(v)d 2 v 

= ~h Lie mW) { I L (v-znv-*>r } dz d * 



and the desired equality, 

(1.8) - / £l = (1 - zw)~\ 

may be verified by expanding the integrand on the left-hand side. 

We point out that the covariance of the limiting resolvent, (1 — zw)~ 2 , 
equals ir x the Bergman kernel for U under that map z — > 1/z (for general 
background on the Bergman kernel, [4] is recommended). The Bergman 
kernel also appears in the covariance of the zero process of a conformally 
invariant Gaussian analytic function on the unit disk, recently studied in [25]. 
In fact, the result of Theorem 1.2 shows that the limiting resolvent <3(z) 
coincides, at least in the sense of finite-dimensional distributions, with the 
Gaussian analytic function 

oo / 1 \ fc+1 

(1.9) G(s) = £vW- 

k=i \ 

in which {Z^} are independent unit complex Gaussians and the prescription 
is valid for \z\ > 1. In this vein, it is also interesting to compare Theorem 1.2 
with the work of [8] on the logarithmic derivative of the characteristic poly- 
nomial of a random unitary matrix. 

Finally, the fact that Theorem 1.2 is stated for \z\ > 4 (and thus the 
analyticity in Theorem 1.1 required for \z\ < 4), whereas the limiting co- 
variance (1.6) and the appraisals (1.7) through (1.9) are on firm ground 
as soon as z lies in the exterior of U, demands explanation. Undoubtedly, 
\z\ > 1 is the correct condition. In fact, for the complex Ginibre ensemble, 
the authors can prove that the corresponding resolvent satisfies Theorem 1.2 
for all \z\ > 1; this work will appear elsewhere. Further, central limit theo- 
rems of the type outlined in Theorem 1.1 are anticipated to hold for test 
functions which are just once continuously differ entiable. The stronger as- 
sumptions made here for more general entry distributions are an artifact of 
the method. Experts on random matrix theory will appreciate the challenges 
in extending techniques from the study of Hermitian random matrices to the 
non-Hermitian setting. 

Since the eigenvalues of M can lie throughout the complex plane, the 
first step in any proof of Theorem 1.2 is to establish an a priori control 
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of the resolvent. With our assumptions granting that the spectral radius 
— |A max (M)| tends to one almost surely, it is natural to attempt to work on 
this event. That is, taking \z\ larger than 1 + e for some e > 0, tr(z — M) _1 is 
well-behaved on the event that |A max (M)| is less than, say, l + e/2. The idea 
would be to prove a central limit theorem for this truncated object, then to 
show that the resolvent times the indicator of {|A max (M)| > 1 + e/2} goes 
to zero in probability in the appropriate norm. However, the mechanics of 
our proof do not allow us to take full advantage of the known estimates of 
the spectral radius. 

The basic idea behind the proof of Theorem 1.2 is to find a martingale 
structure in tr(z — M)~ l by successively conditioning on the sigma-fields 
generated by the columns of M. Given the desired structure, checking that 
the conditions of a standard martingale central limit theorem are satisfied 
requires estimates which compare the resolvent matrix of M and that of 
various rank-one perturbations of M. An example of such a perturbation 
is the matrix formed by setting the fcth column of M to zero, denoted by 
Mfc. The problem which arises is as follows. Having control of |A max (M)| 
does not imply that we have similar control over |A max (Mfe)|. On the other 
hand, the spectral norm of M trivially bounds the spectral norm of Mk for 
all k. Since assumption (1.2.H) implies that the former lies under two with 
overwhelming probability, we can attain uniform control of tv(z — M) _1 , 
tr(z — Mk) , etc., if we are willing to take \z\ larger than four. 

The next section picks up this discussion, describes an overview of the 
proof of Theorem 1.2 and records various estimates needed therein. The 
proof itself then occupies Sections 3 and 4, along with a brief Appendix. 

2. Outline and preliminary facts. Henceforth, we denote the resolvent 
matrix by 

R N ( z )-i = ( z -M)-\ 

As indicated in the Introduction, a characteristic feature of the proof is the 
procedure by which the norm of the eigenvalues of M is controlled. For this, 
we define 

n N = {M : \\M\\ < «}, 

in which || • || is the spectral norm and k is a number larger than two which 
we now fix. Note that for \z\ > 4, there exists a constant K. = K,(k, z) such 
that 

||3jv(z)~ || < £ on the event Ojy. 

Also, from this point, whenever we write \z\ > 4, it should be understood that 
\z\ exceeds four by some fixed e > 0. Thus, the tC in the above bound is also 
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fixed from this point on. We reiterate that ||Af|| < k implies that ||Mfc|| < 
K for any k, with the same implication holding for further perturbations 
of M (in which more than one column is removed, or set to zero). Such 
comparisons constitute the chief reason we are employing the spectral norm, 
rather than the more natural spectral radius, to control the support of the 
empirical eigenvalue distribution. 

Since we essentially wish to work on O^v, it is of course important that 
these events exhaust the space as N — ► oo. In this, we are fortunate that 
the following estimate has been obtained by Geman [13]: under the moment 
assumption (1.2.ii), 

(2.1) P(0$) = o{N- a ) 

for any a > 0. It should be noted that while the proof in [13] is made for real- 
valued random entries, the extension to the complex case is immediate. Also, 
if the entries of M are drawn from a single doubly-infinite array, [31] proves 
that ||M|| — ► 2 with probability one as soon as £7|mn| < oo. However, this 
result does not come equipped with an estimate as sharp as (2.1). 

The stage is now set to describe the proof of the central limit theorem for 
®n(z) =tr3jv(^) _1 — Nz . First, we have the decomposition 

e N {z) = {tvE N {zy 1 l nN -E[tTE N (z)-\n N ]) 

+ (tvE N (z)- l l n c - E[tvE N (z)-\n%]) 

+ (E[trZ N {z)- 1 ]-Nz- 1 ) 

= <5%(z) + & 1 N (z) + <S%(z). 

The Gaussian limit comes out of the first term on the right-hand side, &^(z), 
which is the original resolvent, truncated and recentered. Let Ey. denote the 
conditional expectation with respect to the sigma-field generated by the first 
k columns, mi,m2, • ■ • , rrik of M and write 

N 

(2.2) &%(z) = J2( E k~ E k ^)[trE N ( Z r\n N } 

k=l 

(E Q [-]=E[-]). Clearly, {(^-^^[trS^)- 1 ^^]} form a (bounded) mar- 
tingale difference sequence and we recall the following standard result: 

Lemma 2.1 (Theorem 35.12 of [6]). For each N , let X N:1 ,X Nj2 , ... ,X NjkN 
be a real-valued martingale difference sequence with respect to the filtration 
pN,k> eac h term having bounded second moments. Suppose that for any e > 
and a constant a 2 , 

k N 

(i) lim J2 E l X kk, X N,k\>£}=0, 
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(ii) lim P 



E i X N,k\^N,k] ~ o' 

k=l 



> e 



Then as N — > oo, i/ie distribution ofJ2k=i x N,k converges weakly to a Gaus- 
sian with mean zero and variance a 2 . 

Section 3 proves that appropriate linear combinations of ®^(z) satisfy the 
conditions of Lemma 2.1. Along with a tightness estimate, this will show that 
over a fixed contour C exterior to the disk of radius four, (J5^(z) converges 
weakly to the mean zero Gaussian process with covariance (1 — zw)~ 2 for 
z, w e C. 

The proof of Theorem 1.2 is then reduced to showing that 

(2.3) sup|0k(*)|->O 

z€C 

in probability and that 

(2.4) lim sup\® 2 N (z)\ =0 

for the same given contour C. The proof of (2.4) occupies Section 4; it is 
convenient to dispense with (2.3) immediately. 

Since we have assumed that the densities of the individual entries of M 
are continuous, it follows that sup zgC | tr'E]y(z)~ 1 \ is finite for all iV with 
probability one. (With probability one, the spectrum of M does not intersect 
any predetermined contour.) Therefore, 

pfsupltrS^-^lnc >e) <P(fi&), 
\zec N J 

which goes to zero as N —> oo, by (2.1). To complete the verification of (2.3), 
we must have 

lim supE[\tvZ N (z)~ 1 \,n%]=0. 

This, however, is a consequence of Holder's inequality, (2.1), and the follow- 
ing (rough) estimate, the proof of which is found in the Appendix: 

Lemma 2.2. For any p, 1 < p < 2 and any z EC, 

E\tvE N {z)- 1 \ p <C(p)N 3p / 2+2 . 

The constant C(p) f oo as p j 2, but is independent of z. We note for later 
that the same bound holds if M in Hat(z) -1 = (z — M)~ l is replaced by M^. 

Finally, we conclude this section by recording the following estimate, of 
which we will make repeated use: 
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Lemma 2.3 ([2], Lemma 2.7). For a vector X = (x\, X2, . .., xn) t of in- 
dependent identically distributed (complex) entries x k with E[x{\ = and 
-^[l x i| 2 ] = 1 an d an y X x N (complex) matrix A, it holds that 

E\X*AX - tr A\ p < d {{E\ Xl | 4 tr A*A) p l 2 + E\x x \ 2p tv(A*Af/ 2 ) 

for all p>2 and a constant C\ depending only on p. 

Even more useful will be the related bound: with Ci being a new constant 
still depending only on p, 

(2.5) E\X*AX\ P < C 2 £|xi| 2p ((tr A*A) p l 2 + | tr A\ p ). 

This follows from the above lemma since tv(A* A) p / 2 < (tr A*A) p / 2 for any 
matrix A. 

3. Convergence of the truncated process. We establish the central limit 
theorem for the truncated (and recentered) resolvent (5^. Along with M k 
(or Mke) denoting the matrix formed by setting the kth. column (or both the 
kth and ^th columns) of M to zero, we introduce the additional shorthand 

En^z)" 1 = {z- M k )~\ E NM {z)~ l = {z- Mm)- 1 

(and so forth), for the corresponding resolvent matrices. Likewise, J)jv,fc an d 
ki (and so forth) will denote the events on which M k and M k e have spec- 
tral norms bounded by k. Note that /C still serves as a bound for the spectral 
norm of E^ )k {z)~ l (resp. E^^t{z)~ l ) on the event 0,N,k (resp. Ojv.fe^)- 



3.1. Finite- dimensional distributions. By the Cramer-Wold device, iden- 
tifying the limiting finite-dimensional distributions of 

N N 

<8jU*) = E(^* - E k ^)[trE N (z)-\n N ] = Z N>k (z) 

k=l k=l 

as multivariate Gaussian amounts to verifying the conditions of Lemma 2.1 
for martingale-difference sums of the form 

N ( K ~> N 

M N = J2) ^2 a t z N,k{zi) +/3eZ N>k (z i ) \ = J2 M ^,k 

k=l U=l J k=l 

for any fixed K, any choice of z^s on C and complex (o^,/^) such that the 
resulting object is real. As a first step, A4n is simplified through a series of 
reductions in each of its terms. 
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Lemma 3.1. As N — > oo, Mn either converges in distribution or fails 
to along with a related martingale- difference sum Mn- The latter is defined 
by replacing each appearance of Z^ k {z) i n -M-N,k with 

ZN,k(z) = E k \e k E NM (z)~ m k ,Q, Ntk ] = e k E k [E Ntk (z) ,£l N ,k\ m k, 

where e k is the canonical unit vector in M. N . In other words, the difference 
of Mn and Mn goes to zero in probability. 



The point of Lemma 3.1 is the introduction of independence in each term, 
particularly that of the vector m k and the matrix E k [E^ >k (z)~ 2 , fijv^]. With 
that established, the desired conclusion comes down to the following: 



Lemma 3.2. The martingale- difference sequence Mn satisfies the con- 
ditions of Lemma 2.1. In particular, the limiting covariance (1.6) in Theo- 
rem 1.2 is seen from the fact that 

N 

(3.i) E^-it^W-* E ^(l-zizjr 2 

k=l l<i,j<K 
in probability as N — > oo. 



Proof of Lemma 3.1. Since (E k - E k ^i)[tvE N:k (z) 1 ] = 0, 

N N 

E Z N,k( z ) = ^2( E k ~ E k ^ 1 )[tliE N (z)~ l - Sjv,fc(2;) _1 ),Ojv] 
k=l k=l 

(3 - 2) 

-J2( E k~ E k ^)[tTE N)k (z)-\n%]. 

k=l 

By Holder's inequality, estimate (2.1) of P(f2^) and Lemma 2.2, we may 
conclude that 

in which the constant Co is independent of the index k. That is, the second 
term on the right-hand side of (3.2) tends to zero in probability as N — > oo 
and thus may be neglected in questions of distributional convergence. For 
the first term, we bring in the usual resolvent identity along with the basic 
fact that 

(3-3) (A + ve T k )- l v= A ~^ V 

(l + e^A L v) 
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(3.4) 



for any invertible A and vector v. Applied in succession, this yields 
tr(E N (z)~ 1 - En^z)' 1 ) = tr(E N (z)~ 1 (m k el)E N)k (z)~ 1 ) 

= el(E Njk (z)~ 1 E'^(z))'m k 

<%E N , k (z)-*m k 
^ + elE Ntk (z)- 1 m k 

= Sk(z)e k E N ^ k (z)~ 2 m k . 

Next, presuming control on the norm of E^^ k {z)~ l and the fact that m k is 
mean-zero and normalized as in £7||mfc||| = Yli=i E\rri£ k \ 2 = 1, it is antici- 
pated that e^EN^ k (z)~ l m k tends to zero as N — > oo. That is, we wish to 
replace 5 k (z) with one in each appearance of the sum (3.2), and conclude 
that 



N 



(3.5) ^2(E k - E k -i)[e k E Ntk (z) 2 m k - S k (z)e k E N>k (z) 2 m k ,Vt 
k=i 



N 







in probability. Expanding a little, we find 

1 - S k (z) = -e k E N}k {z)~ 1 m k + (e k E NM (z)~ l m k ) 2 

- (e k E N ^ k (z)~ l m k f + 5 k (z)(elE N)k (z)~ 1 m k )' i . 

The sum in question has been broken into four, the first three terms of which 
may be bounded in mean-square as follows. With p = 1, 2 or 3, 



E 



N 



Y^( E k ~ E k -i)[{e k E N:k (z) 1 m k ) p (e k E N:k (z) 2 m k ),n N ] 



k=l 



N 



(3.6) 



< E[\e k E Njk (z) 1 m k \ 2p \e k E Ntk (z) 2 m k \ 2 ,n Ntk ] 



k=l 
N 



<AY,{^rntA N)k (z)m k \ 2p E\m* k B Nyk (z)m k \ 2 }^ 2 , 
k=l 

in which Aj\r^ k (z) and Bjv,fc(.z) are the Hermitian matrices 
A-N,k( z ) = (H A r ifc (z)" 1 )*(efcefc)HAr ifc (z)" 1 lo Jvfc 

and 

&N,k( z ) = {^N,k{z)~ 2 )*{e k el)E N , k {z)- 2 lu. Nk . 

Having replaced £In with Ojv^ D in line two of (3.6), note that m k is 
independent of both Ajv k (z) and Htf k (z). This makes way for an application 
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of Lemma 2.3 [actually, its consequence, (2.5)]. We have 
E\rn* k A N>k (z)rn k \ 2p 

<dA^£[(trA^(z)n 
(3.7) =C 1 N" 2 PE[(el(E N , k (zr\E^ k (zy 1 r)e k r,n N , k } 

<C 1 N- 2p E[\\E Ntk (zr 1 f P ,^N,k] 
< dlC 2p N- 2p , 

where the finiteness of the constant C\ relies on the assumption that the en- 
tries of M have bounded moments. Since the expectation involving B^/t^) 
in (3.6) may clearly be estimated in the same way, it follows that the last 
line of that display is bounded above by a constant multiple of N~ p . 

As for the last sum resulting from the expansion of (3.5), notice that on 
the set Qn, the simple bound 

\5 k (z)elZ N:k {z)- 2 m k \ < | tr H 7V (^)" 1 | + | tr H^z)" 1 ] < 2NK 

may be read off directly from the identity (3.4). It follows that 



E 



N 



^2(E k - E k ^ 1 )[(elE N)k (z) 1 m k ) A 5 k (z)elE Ntk (z) 2 m fc ,0 



IV. 



k=i 



N 



(3.8) 



< 16iV 2 /C 2 J2 E\m* k B N>k (z)m k \ 



k=l 



This completes the verification of (3.5). 

To produce the aforementioned form of Mjy,fc requires two more steps. 
First, we have 



E 



N 



^2(E k - E k -i)[e k E N 2 k (z)m k ,n N>k ] 

k=l 

N 

~J2( Ek ~ E k -i)[elE N 2 k (z)m k ,n N 

N 

^2(E k - E k _ 1 )[el'B.] f 2 tk (z)m k ,n Njk n 



k=l 



E 



k=l 



<2Y J {E[\elE N)k {zr 2 m k \\n Ntk ]) l, \P{^ c N )f li 

k=l 
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N 

= 2Y J {E\mlB N ^z)m k \ 2 ) l l\P{n c N )fl' 
k=l 

which goes to zero faster than any polynomial power of N [recall (2.1)]. It 
remains to observe that 

E k _ l [el'E N)k (z)~ 1 m k 9. Ntk ] = e^ fe _i[H A r ifc (z)~ 1 , 9. N>k }E k _ x [m k ] = 0, 

the independence being the point of replacing Qjv with £l]y,k termwise. The 
proof is thus complete. □ 

Proof of Lemma 3.2. Verifying the first condition of the martingale 
central limit theorem (Lemma 2.1) is trivial. In the previous proof, we 
have already made repeated use of the fact that E\Zn k (z)\ is less than 
2E\m* k ~f$N^ k (z)m k \ 2 = 0(N~ 2 ) for arbitrary z €C and all k. Hence, 

JV N 

£ E[M 2 Ntk , \M N , k \ >e]< e~ 2 £ E[M%, k ] 

k=l k=l 

N K 

< c 2 ci E E(N + mfE[\z N ^z l )\ i \ 

k=ll=l 

tends to zero as N — > oo. 

Next, by the linearity of M.N,k in the Zjvfe's, the covariance formula (3.1) 
follows from checking that 

JV 

(3.9) E k ^[Z N , k {z)Z N , k {w)] = 
k=i 

and 

N „ ^ l 

(3.10) ^^fe^H]^—^ 
fc=i ^ z ^ 

in probability as iV —> oo. For the equality, using the independence of m k 
and the feth row of 

&N,k{ z ) = E k [ZN )k (z)~ 2 , fijv,fc]> 

we find that 

TV 

^ Sfc~l[^iV,fc(2)^iV,A : (^)] 
k=l 
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N 

= E k ^i[(elB N)k (z)m k )(elB N:k (w)m k )] 
k=l 
N 

= X] [SN,k(z)}ki[SN,k(w)}kjE[mi k m jk ] = 0, 

k=ll<i,j<N 

as both E[rriij] and E[mfj] are assumed to vanish. On the other hand, with 
conjugates present in the second half of (3.10), the computation becomes 

N 



Y E k _ x [Z Njk (z)Z N>k (w) } 

k=l 

1 N 



N k=i 



N 



= -^lt e k( E k[(z- M k y 2 ,n Ntk ]E k [(w - M fc *)- 2 ,0 7V ,fc])e fc , 
k=i 

in which we have temporarily removed the notation and written everything 
out in full for the sake of clarity. 

The next step is to replace the question of the convergence of (3.11) with 
that of its antiderivative. In particular, what we proceed to prove is that 

1 N 

JjY e I( E k[^N,k( z ) 1 ^N,k] E k[{-N,k( w ) l T^N,k]) e k 
k=l 



(3.12) 



log( 

zw 



in probability. That this is sufficient follows from a theorem of Vitali (see, 
e.g., [29]): If {.fjv(z)} is a sequence of uniformly bounded analytic functions 
on a domain B C C, converging at a collection {z} containing a limit point 
in B, then there exists an analytic f(z) such that /n(z) — ► f(z) and f' N (z) — > 
f'(z) throughout B. Denote the left-hand side of (3.12) by 4>n{z,w). This 
function is analytic in z {\z\ > 4) for any fixed w (\w\ > 4); in the same region, 
it is analytic in w for any fixed z. That 4>n{z-,w) is bounded independently 
of iV follows from the fact that 

(3.13) ll^tH^^)- 1 ,^]!! 2 ^^^!!^^^- 1 !! 2 ,^]^^ 2 . 

Now, if the convergence (3.12) takes place for given z and w, we can fix 
a collection {z k ,w k } with limit point in {|z|,|u)| > 4} and a subsequence 
N' — > oo over which the convergence is almost sure, simultaneously, for all 
points (z k ,w). On such a realization, the theorem of Vitali can be applied 
twice, one coordinate at a time, to obtain the convergence of 4>N an d its 
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mixed partials over the full domain. Since the conclusion is independent 
of both the subsequence and the particular realization, (3.12) does indeed 
imply the latter half of (3.9). 

Moving on, we bring in some additional shorthand, 

A N ,k(z) = E N)k (z)~ 1 ln N-k , A N ,ke( z ) = s Af,M(^)" 1 ln JV , w , and so on, 
and focus on a fixed term of (3.11), 

(3.14) T Ntk (z,w) = elE k [A N)k (z)}E k [A N , k (w)*]e k . 

The method is rather laborious. First, successive columns of M are removed 
[effectively replacing An^{z) with An,kt{ z )i etc.] in order to invoke the ba- 
sic estimate of Lemma 2.3. Afterward, the removed columns are effectively 
restored, arriving at an equation for T]\r jk (z,w) which is approximately de- 
terministic for large N. 

Invoking the matrix identity 



(3.15) Z N)k (z) 1 = -I +-^2E N)k (z) 1 m e ej 

z z i^k 



T 



(3.14) is expanded: with P k denoting the conditional probability with respect 
to mi, ... , mk , we have 

T N , k (z,w) 

= —Pi(n N>k ) + — Y^(e k T E k [A N , k {z)me\ei T e k )P k (Sl 



(3.16) 



. N,k) 

ZW ZW — 

+ — y2e k e£E k [rn}AN,k(w)*]ek 

+ — e[ E k [A N ,k{z)me]ei T ejEklmj* A N ,k(w)*]e k . 

ZW t^kj^k 

This, in turn, may immediately be reduced to 

T N ,k{z,w) = — + el E k [AN,k(z)rn£]E k [m e *A N ,k(w)*]e k 

zw zw f^ k 

(3.17) 

+ S Nfk (z,w). 

Here, S^ k (z,w) = (zw)~ l {l — P^(J7jv,ifc))) which is easily seen to go to zero 
in L . Note also that the second and third terms of (3.16) are identically zero: 
both contain e^e^ = in each £ ^ k term of their defining sums. The reduc- 
tion of the last term of (3.16) to the sum present in (3.17) follows similarly. 
As for that sum, we wish to replace the -4tv,A;(')' s with the corresponding 
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•A.N,ki{-)- The procedure and goal, which is to produce independence between 
the resolvent matrices and the columns of M, follows the outline of the proof 
of Lemma 3.1. 

First, note that for £ ^ k, 

'E N ^{z)~ 1 m(, = 5u{z)E Nj u{z)~ l m^ 

in which 

hi{z) = (1 + ^ T '^ NM {z)~ l m t y 1 = (1 + e M (z)y l 
[recall (3.3)]. Now write 

^2elE k [A Ntk (z)m e ]E k [m e *A Ntk (w)*}e k 

i^k 

= ^2 e k E k[^N,ke(. z )~ lm e^N,k\ 

(3.18) 

x E k [mi*(E NM (w) )*,fijv,fc]e fc + 

= elE k [A N ,k£( z )}( m i m i*) E k[-^N,k£( w y] e k + £j\r,fc + £ N,ki 
Kk 

where E k \AN, k £{-)mg\ = for £ > k is used in the last line. Sj^ k = £jj k (z, w) 
and £ff k = Sfj k (z, w) are the sums over £ ^ k of 

£ N,k,A z ^) = E k [{6 ki (z) - l)efcS A r,^(^)~ 1 m £ ,0 7 v,fc] 
x E k [mi '(En^w)' 1 )* 'e k ,Q Njk \ 
+ E k [e k "'~ NM (z)~ l rn t ,Q, Ntk ] 

x E k [m(*{E NM {w)~ l )*e k (5 M {w) - l),^N,k] 



(3.19) 



+ E k [(5 k e(z) - l)e k E N)H (z) 1 me,Q N ,k] 
x E k [m e * '(En^w)' 1 )* e k {5 M (w) - l)Q N ,k] 

and 

£ N,k,e( z ^ = -E k [elA NM (z)?n e ]E k [m e *A NM (w)*e k ,n% !k \ 

(3.20) - E k [elA NM (z)m e ,n% jk ]E k [m e *A NM (w)*e k ] 

+ E k [elA NM (z)m e ,n% jk }E k [m e *A Niki (w)*e k ,n% )k \. 
Similarly to (3.7), that is, by Lemma 2.3, we have 

^|^fc[efcSjv,fe«(-)~ lm ^. n JV,fc]| 2 < E[\elE NM (-)~ l m,,\ 2 ,VL NM ) 

(3.21) = E\m}A N ^{-)*e k el A Njk> e(-)m £ \ 

< c^- 1 
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for • = z,w. Next, note the equality 

(S M (-) - l)(e^'r,]\f^e{')~ lm e) = (-£«(•) + £M{-) 2 5 k i{-)){elE NM (-)~ l m t ) 

= -{ee T E NM (-)~ 1 m e )(elE NM (-)~ 1 m e ) 
+ (ei T E N)k i(-)~ 1 mi) 2 (elE N:k (-)~ 1 m i ), 

the identity 5 k i(-)EN : ke(')~ lm e = ^N,k(')~ lm e having been reinvoked in the 
last step. It then follows by the same type of reasoning used in (3.21) that 

(3.22) E\E k [(S kt (-) - l)e^ ifc ,(-)~V,fW-]| 2 < 2Xi + 2T 2 , 
where 

(3.23) Xi = J B[|e^ T 3 A r,H(-) _1 ^| 2 |efc5 A r,w(-) _1 ^| 2 ,^iV ) fe] 

< E[\e£ T ANM(') m A 2 \ e k^NM(') m i\ 2 ] 

< C 2 N~ 2 , 

after an application of the Schwarz inequality and 

(3.24) l 2 = E[\e l T E NM {-r 1 m,\ i \elE N ^ 1 m l \ 2 ^ Ntk ] 

<JC 2 E[\e e T A NM (->M 4 \\m e \\ 2 2 } 

<IC 2 (E[\\me\\i}) l / 2 E[\e e T A N M>M\ 8 } 1/2 
< C 3 N~ 2 . 

Here, || • ||2 denotes the usual i -norm. 
Combining (3.21) and (3.22)-(3.24) gives 

E\£ 1 N)k (z,w)\=0(N- 1 / 2 ). 

On the other hand, (3.21) plus the sharp decay P(£l% k ) is enough to produce 
a similar bound on the L 1 -norm of £% k (z,w). 

Returning to (3.18), since E[m£m}] = (1/N)I, it is natural to claim that 



(3.25) 



T7 elE k [A NM {z)]E k [A NM (w)*]e k + £%^ k {z, w) 



N Kk 



with a small error; we show below that 

(3.26) E\£l k (z,w)\=o(l). 
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At this point, the endgame becomes clear. If we could now restore the last 
column removed from the resolvent through an estimate along the lines of 

E\el(E k [A NM (z)]E k [A NM (w)*} 

(3.27) 

- E k [A Njk (z)]E k [A Njk (w)*})e k \ = o(l), 

the combined outcome of this statement with (3.17), (3.18) and (3.25) would 
be that 



T N , k (z,w) 



1 



+ 

zw zw 



1 fk-l, 



N 



T N , k {z,w) +£ N>k (z,w), 



in which 



E\£ N)k (z,w)\ = o(l), 



uniformly in k, 1 < k < N. This is to say that the object under investigation, 
the left-hand side of (3.12), satisfies 



N 



z,w) 



k=l 



1 N 

-Y 



zw 



k-1 

N 



< 



\zw\ 



zw 



1 



N 



w 



k=l 



and so tends to zero in probability. The advertised limit follows from iden- 
tifying the Riemann sum. 

Returning to the bounds (3.26) and (3.27), we begin with the second. The 
arguments used in controlling the error £jy , (z,w) above will show why 

E\elE k [A N ,ki( z )} E k[-AN,ki(. w Y} e k 

~ e k Eki^NMi 2 ) -1 ' ^■N,k]Ek[(^N,ke( w )~ 1 )* ^ ^N,k]ek\ 

tends to zero faster than, say, iV^ 1 . Next, to restore Hjv,fc^(-) _1 to H7v,fc(") _1 > 
note, for example, that 

E\E k [el (3jv,jfc(^) _1 - E N , k e(z)~ 1 ),tt N , k ]E k [(Z N , k (w)~ 1 )*e k ,Q Ntk ]\ 2 

<E[E k [\elE NM {z)~ l m^{E N ^ k {z)- l fe4tn N ^} 

x E k {\\{E NM {w)- l )*e k \\l,n NM }} 

< JC 4 E\elAN,ke( z ) m e\ 2 
= 0(N- r ). 

It will follow that the right-hand side of (3.27) is of order N~ 1 ' 2 for any k,£. 

Finally, to verify (3.26), we take the mean-square of that error term, 
written as 



N 2 



E 



^{mlE k [ANM( w )*]( e k e l) E k[ANM( z )]^e- 

Kk 
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(3.28) 



(3.29) 



elE k [A N ,ke( z )] E k[^N,ke( w )*} e k) 



^(m^D N)k i(z,w)me - tr(D NM (z,w))) 

Kk 



where rhp = v Nrri£ are the columns of the unnormalized matrix. Again, by 
Lemma 2.3, we have that 



(3.30) 



E\m* e D Ntki (z,w)rhe - tT(D N ^i(z,w))\ < C 4 



for a fixed constant C4. Therefore, the contribution from the diagonal terms 
of (3.28) is of order iV _1 . Consideration of the off-diagonal terms is compli- 
cated by the dependence between E k [AN,ki(')} and E k [Ai\[ jk j(-)], forcing us 
to remove yet another column to arrive at the required estimate. Note that 
with j, I < k, the typical cross term satisfies 

E[(rh* £ D NM (z, w)m e - tv(D NM (z, w))) 
x (m*D Nikj (z,w)fhj - tv(D Ntkj (z,w)))] 
= E[{fnlD NM (z, w)rh(, - tv(D NM (z, w))) 
x (rh*(D Ntkj (z,w) - D N)kj t(z,w))rhj 
- tr(D N)kj (z,w) - D N:kji (z,w)))}, 

where in Dj^ >k ji(-, •), each Azv,/y(-) is replaced by AN,kji(')- After two ap- 
plications of the Schwarz inequality and (3.30), the required estimate will 
follow if 



(3.31) 

and 

(3.32) 



E\fh*(D N>kj (z,w) - D Njkje (z,w))mj\ 2 = o(l) 



E\tr D Ntk j(z,w) - tiD N:kj£ (z,w)\ 2 = o(l) 



for iV — > 00. However, Lemma 2.3 bounds the left-hand side of (3.31) in terms 
of the left-hand side of (3.32). It remains to observe that after unraveling the 
notation, (3.32) is really the same as (3.27). We need to compare combina- 
tions of the resolvent with two columns removed, ejT E k [AN,ki(z)]E k {AN,k£(w)*]e k , 
to the same object in which a third column has been removed. (3.27) com- 
pares that same combination with one column removed to its image after 
a second column has been removed. An identical procedure used to bound 
(3.27) will then show that (3.32) is 0{N~ 1 / 2 ). The proof is thus complete. 
□ 
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3.2. Tightness estimate. The previous section established the conver- 
gence in distribution of any finite collection {®^(zi), (5^(22), ©jv(<Zfc)} 
of the truncated and recentered resolvent 

&%(■) =trE N (-r 1 ln N - E[tiE N (-)-\n N ] 

TV 

= (E k - Ek-i) [tr E N (z)" 1 , n N ] , 
k=i 

sampled at any points Z£, 1 < £ < k, exterior to the disk of radius four. The 
standard method for extending such finite-dimensional convergence to the 
convergence of <&%(z) in the space of continuous processes over a prede- 
termined contour C is to check the Arzela-Ascoli criteria (see, e.g., [5]): 
verifying the estimate 

E\®° N (z)-<5%(w)\ a <C\z-w\ 1+ P 

for any a, f3 > and all z, w £ C is enough. 

In the present case, it is convenient to recenter once more before preceding. 
Note that 

TV 

]T(£ fc - £ fc _!)[tr -^(z)" 1 ] = 
fc=i 

for any z, while for any C exterior to the disk of radius four, 

TV 

supJ2( E k - E k ^)[tiE Ntk (z)-\n%]^0 
z&c k=i 

in probability. The proof of the latter is the same as that for sup 2gC |0]^(z)| 
— ► in Section 2. It follows that the tightness of &%(■) in the space of 
continuous functions on a given contour C is equivalent to that of 

TV 

&n(z) = T,(Ek- tffc-OIMHjvCz)- 1 - Sjv.kCz)- 1 ),^]. 

k=l 

The following estimate, along with Lemma 4.1 proved in the next section, 
completes the proof of Theorem 1.2: 



Lemma 3.3. It holds that 



E 



®%(z) - <5%(z) 



z — w 



2 



with a constant C independent of both N and any choice of z, w such that 
\z\, \w\ > 4. 
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Proof. It is no surprise that considerations are much the same as those 
behind the proofs of Lemmas 3.1 and 3.2. We at once have that 

N 
k=l 

x[tr(E N (z)~ 1 E N (w)~ 1 - Ejv jfc (2) _1 Ejv )fc (iu) -1 ),fijv]. 
Now, again with 6 k (-) = (1 + efcSjv,fc(-) _1?T &fc) _1 j we write 

E N (z)~ 1 E N (w)~ 1 - E^zy^N^wy 1 

= (E^zy 1 - E^y'XENiwy 1 - e^) -1 ) 

+ (E N (z)~ 1 - E Ntk (z)~ 1 )E Ntk (wy 1 

+ E N:k (zy 1 (E N (wy 1 - Eu^wy 1 ) 

= S k {z)5 k (w)(E N>k (zy 1 (m k e k ) 

x E Ntk (zy 1 E N:k (wy 1 (m k el)E N (wy 1 ) 
+ 5 k (z)(E N:k (zy 1 (m k el)E Ntk (zy 1 E N>k (wy 1 ) 
+ Sk{w)(E N>k (zy 1 E N:k (wy 1 (m fc efe)Hjv,jfe(^) _1 ) 

and so 

tr(Hiv(^) _1 Ejv(i(;) _1 - Ejv,fc(.z) _1 E;v,fc(w) _1 ) 
= S k (z)5 k (w)(elE Njk (zy 1 E k (wy 1 m k ) 2 

(3.33) 

+ 5 k (z)(elE N)k (zy 2 E N>k {wy l m k ) 
+ 5 k (w)(e k E N>k (wy 2 E N ^(zy 1 m k ). 
We therefore wish to show that 

(3.34) ElS^Skiw^elEN^zy^kiw)-^) 2 , n N \ 2 < CiiV" 1 
and 

(3.35) EMz^elE^zy^kiwyW),^ 2 < C^ 1 

for arbitrary k, z and w [the last two terms of (3.33) are the same up to 
exchanging z and w]. 

We provide the full details behind (3.34) only, this being the more elabo- 
rate and both assertions requiring the same strategy. Further, the strategy 
is familiar. In particular, the reader will recognize that the desired bounds 
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follow directly from Lemma 2.3 if 6k(z) and 5 k (w) can be replaced with 1 
inside both expectations: 

E\{el'E N ^(z)~ 1 'E N ^{w)~ 1 m k ) 2 ,Vt N \ 2 

< E[\ml(E NM (z)~ l E N ^ k (w)~ l )*{e k el){E NM (z)~ l E N ^ k (w)~ l )m k \ 2 , 

^N,k] 

(3.36) 

< C 2 N~ 2 E[\e T k (H JVifc (z)- 1 H JVifc H- 1 )(H^ fc (z)- 1 H^H- 1 )*e fc | 2 , 

^N,k] 

< C 2 JC 8 N~ 2 
and, likewise, 

E\elE N)k (z)~ 2 E Ntk (w)~ l m k ,tl N \ 2 

< ^[|mfc(H 7 v,fc(z)" 2 H 7V ,fc('w)~ 1 )*(eA : efc)(H A r ifc (z)~ 2 H A r ifc (w)~ 1 )m fe |, 

Recall the notation 5 k (z) = (1 + e k (z))~ l and set 
5 k (z)5 k (w) = (1 + e k (z) + e k (w) + e k {z)e k {w)y l = (1 + -y k (z, w))~ l . 
The key step in proving (3.34) is to notice the a priori bound: on the event 

\5 k (z)5 k (w)(elE N ^ k {z)~ 1 E N ^ k (w)~ 1 m k ) 2 \ 

= Itr^z)- 1 - Hyv.fc^r 1 )^^)" 1 - H^fcH" 1 )! 

(3.37) 

< ivd!^)- 1 !! + HH^^-^Ddi^H- 1 !! + p^H- 1 !!) 

< AK 2 N. 

Next, we know by now that E[\e k (-)\ 2p , f2jv,fe] < C^N~ P for p > 1 and thus 
there is a constant C5 such that 

(3.38) ^(z,^,^] <C 5 N-f. 

If we now write 

5 k (z)5 k (w) = 1 - 7fe(z,w) +7fc(z,if) 2 - 7/ c (z,'u;) 3 (5fe(2;)(5fc(w) 

in the left-hand side of (3.34), successive applications of the Schwarz in- 
equality, along with (3.36) and (3.38), will show that each of the resulting 
first three terms are of order N~ 2 . For the fourth term, the bound (3.37), 
followed by the appraisal E[\^y k (z, w)\ 6 , Cl^] = 0(N~ 3 ), will produce the or- 
der A r_1 in (3.34). For (3.35), the argument is as follows. Playing a similar 
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role to (3.37), the fact that 6k{z)(Epf ! k(z)~ 2 Ek(w)" 1 rrik)ln N is less than a 
multiple of N is noted. Then 6f.(z) is itself expanded to third order with 
each term being handled in an identical manner to that just described for 
(3.35). The proof is thus complete. □ 

4. Equivalence of the means. We prove the following: 

Lemma 4.1. It holds that 

lim ^[tr-Tv^)" 1 ] - Nz^l = 0, 

uniformly for z belonging to any contour exterior to the unit disk. 

In the notation of Section 2, this reads as — ► 0. Therefore, the 
central limit theorem is valid for linear statistics centered by the asymptotic 
mean as well as the iVth level mean. Note that the above result requires 
only that \z\ > 1. As we have discussed, this is certainly the appropriate 
(i.e., optimal) condition for Theorem 1.2. While the technique used to prove 
the CLT rests on comparing rank one perturbations of the random matrix 
M, the present computation requires no such comparisons. Thus, rather 
than relying on bounds on the spectral norm to provide a priori control, we 
are able to use the spectral radius which is more natural to the problem. 



Proof of Lemma 4.1. Denote 

A 7v = max{|A| : A is an eigenvalue of M}. 
Along with (2.1), Geman [13] has also proved that for any p > 1 and a > 0, 

(4.1) P(A N>P )=o(N~ a ). 

(Once again, [13] discusses the case that the entries of the matrix are real, 
but the extension to complex entries is effortless.) Also from [13] [see the 
display after equation (6) of that paper] , we record the bound 

(4.2) E[A p N ]<CtN 2 p 

for a constant C\ and p allowed to grow with N, up to p = 0(\nN). 

Now, fix a positive p larger than one but less than \z\. The bound (4.1), 
together with the often-used Lemma 2.2, implies that 

E[trE N (z)- 1 ,A N >p]^0 

as N — ► oo. On the opposite event, E^(z)~ 1 may be expanded inside the 
expectation as in 

E[E N (z)-\A N <p] 

= E -m E ^ Mk > An ^ p] + z^ e M*n(*T 1mp )^n < p\- 

k=0 z z 
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The first term in the sum on the right-hand side reads z~ l {N — NP{\n < 
p)) = Nz~ x + o(l), in which we see the desired result. Next, note that the 
remainder term is bounded by N(\z\ — p)~ 1 (p/\z\) p and thus goes to zero 
for p = p(N) = C^lniV with a sufficiently large constant C2. Further, on 
account of (4.2), 

E[trM k ,A N >p] < NE[A% } 1/2 P(A N > p) 1/2 < CsiV" 1 

for all N sufficiently large and all k < (C2/2)lniV. Therefore, it remains to 
prove an estimate of the type 

(4.3) lim V E[tvM s ] = 0. 

l<s<CilnJV 

By the assumptions that E[rrike] = ^[m 2 .^ = 0, we have 

\EtxM s \ < y E\m ei £ 2 me 2 e 3 ■ ■■m ig i 1 \, 

in which J2* denotes the sum over those choices of indices {£±,£2, ■ ■ ■ ,£ s } 
such that any index pair (£j,£j + i), or entry type m^^ +1 , which appears 
must do so at least three times. This object may be bounded in a crude but 
sufficient manner as follows: 

E\m ei e 2 me 2 i 3 ■ ■■m isll \ 



IV3J 



fc=l niH h"fe=s 



E -^-^ :N k (E\m u rE\m u r---E\mn 

jL — i £ — i "i Tin •■•Bl 



The outermost sum indicates that k = 1, . . . , [s/3\ distinct entry types can 
be present in any string being summed over. The second sum is over all 
possible numbers of each entry type: n\ type 1 entries, ri2 type 2 entries, and 
so on. We allow for all nonnegative, ordered fc-tuples ni, ri2, ■ ■ ■ , n& summing 
to s, which results in an overestimate. The next factor, s\/n\ \ ■ ■ ■ n^!, assigns 
each entry from left to right one of the k types. Afterward, each type must 
be labeled by choosing the corresponding pair of indices. That is, the first 
entry in the string might have been assigned "type 3" at this point, but type 
3 can be m^2, or 7717,4, and so on. Note that not all such unrestricted type 
assignments will support a consistent labeling of indices; we have another 
overestimate. As for the number of ways to label the indices, the factor 
N k is an obvious upper bound. Moving from left to right, each time a new 
(unassigned) entry type is encountered, one of N possible choices for its first 
index is chosen. After the first index corresponding to each particular type is 
chosen, all of the second indices are determined by their neighboring entries. 
Again, we allow for the overcount produced by a number of inconsistent 
labelings which result from this procedure. The final product of expectations 
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requires no explanation and by our assumptions [see (1.2.iii)], it is bounded 
as in 

k 

E\mu\ ni E\ mil \ n2 ■ ■■E\m 11 \ Uh < N~ s / 2 JJ(n < ) an < < N~ s / 2 s as 

e=i 

for some a > 0. It follows that 

[s/3] 

\EtrM s \ < N~ s / 2 s as N k k' < N- s / 6 s {{a+1)s+1) . 

k=l 

This clearly goes to zero as N — > oo with s = OQnN), and it does so fast 
enough that the sum over all such s will also vanish [recall (4.3)]. The proof 
is thus complete. □ 

APPENDIX 

We include here a proof of Lemma 2.2. Important use is made of the 
following fact which was established in [1] for a similar purpose: 

Proposition A.l (Corollary A. 2 of [1]). Let the vector X = (xi,x 2 , ...,x ri 
consist of independent complex entries such that the joint density of the real 
and imaginary parts of each xg are uniformly bounded by a constant C. Let 
ai, a.2, ■ ■ ■ , CKfc be k orthonormal complex unit vectors. Then the joint density 
of the complex numbers X*a±,X*a2, ■ ■ ■ ,X*a k is bounded by C 2k n 2k . 

Proof of Lemma 2.2. For any N x N matrix A, \ J2k=i ^k(A)\ is dom- 
inated by EfcLiP*^)fcfc] 1/2 and so 

N 
k=l 

when p > 1. The notation on the right indicates (the p/2 power of) the 
kk entry of the inverse matrix of En(z)*Ej\[(z). Next, let the kth column 
of Hjv(-z) be denoted by £k( z )- Also, let E N ^(z) denote the iV x (A — 1) 
matrix formed by removing £k( z ) from 3n(z). Then 

(A-i\ (~ M*n r^-i - det ( E N(k)(z)*Z N{k) {z)) 1 
(A.l) {a N {z) ^N{z)) kk - ~ F77^77T77> ' 

aet{r. N {z) z<n{z)) gb(A) QkZk\z) 

where 

Qk = Qk(z) = L - E {k) (z)(E^ ) (z)*E (k) (z)y 1 E^ ) (z)* 

is the projection onto the orthogonal complement of the space spanned by 
columns 1, 2, . . . , k — 1, k + 1, . . . , N of Ejv(z). The second equality in (A.l) 
rests on the fact that det(E]y(z)*EN(z)) = det(E N ^(z)*E N ^ k )(z)) times 
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{£,k(z)* : Qk£k{ z )) > which is an instance of the Shur complement formula (see [17], 
page 21). 

Next, since the entries of M, and so Xn(z), are continuous random vari- 
ables, each Qk is rank one with probability one. Hence, there are complex 
unit vectors 7^ such that = Jklk- Summarizing, we have the bound 

N r 1 



£7|trS JV (z)~ 1 | p < N^Y, E 



k=l 



Mz)*j k \ p 



Note that the vectors £k(z) and 7^ are independent. By assumption, we also 
know that joint densities of the real and imaginary parts of each entry of 
the scaled vector y/N£k(z) are uniformly bounded. From Proposition A.l, it 
follows that the density of the complex random variable (\fN£,k( z ))*lk con- 
ditional on 7^ is bounded by a constant multiple of iV 2 . Since that constant 
is independent of the index k and z, we conclude that 

< C{p)N 3p / 2+2 , 
which completes the proof. □ 
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